Stress induced dislocation roughening — phase transition in Id at finite temperature 
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We present an example of a generically forbidden phase transition in Id at finite temperature - 
stress induced and thermally assisted roughening of a superclimbing dislocation in a Peierls potential. 
We also argue that such roughening is behind the strong suppression of the superflow through solid 
*He in a narrow temperature range recently observed by Ray and Hallock (Phys.Rev. Lett. 105, 
145301 (2010)). 
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PACS numbers: 67.80.bd, 67.80.dj, 67.80.-s 

Strong interest in the supersolid state of matter in free 
space [l| has been revived by the recent discovery of the 
torsional oscillator (TO) anomaly in solid ^He [2[. While 
finding no supersolidity in the ideal ^He crystal, ab ini- 
tio quantum Monte Carlo simulations did find that some 
grain boundaries Q , dislocations [1, Q or crystal bound- 
aries support low-d superfluidity spatially modulated 
by the surrounding lattice. In principle, a percolating 
network of superfluid dislocations could explain the 
TO anomaly if the dislocation density is 3-4 orders of 
magnitude higher than it is expected to be in a slowly 
grown and well annealed crystal. Consistent with such 
expectation is also a very small rate of the critical super- 
flow through solid ^He (occuring presumably along dis- 
locations with superfluid cores) observed in the UMASS- 
Sandwich experiments [1, Q . Thus the nature of the TO 
anomaly in solid '^He remains unclear. 

In the present work we focus on the very unexpected 
feature of the UMASS-Sandwich experiment - the 
strong suppression of the supercritical flow rate Vcr (by 
about 3-4 times!) and then its recovery in a narrow range 
of temperatures. Such a feature occurs well below (about 
10 times) the flow onset temperature Tq « 0.5 — 0.6K 
Here we are proposing an explanation within the model 
of superclimbing dislocation [5|, biased by an externally 
imposed chemical potential fi which generates a stress on 
the dislocation core and, thus creates spontaneous jog- 
antijog pairs. 

Jog-antijog pairs as quantum objects can be created 
spontaneously by a macroscopically small stress a > ac (x 
1 /L applied to a superclimbing dislocation of length L - 
analogous to the creation of kink-antikink pairs along a 
stressed gliding dislocation [l^, [HI . We have found that 
such an instability leads to a flrst-order phase transition 
even at finite temperature T between two phases of the 
dislocation - smooth and rough. This transition is in an 
apparent violation of Landau's argument " no phase tran- 
sitions in Id at finite T" [l^]. However, we argue that 
the locality of order parameter(s), which is essential for 
the validity of Landau's argument, is not present here. 
Consequently, in a sharp contrast with conventional Id 
systems, where any macro scop ic characteristic length de- 
creases with increasing T [ij] , a typical scale Lh for the 
onset of the hysteretic behavior of the superclimbing dis- 



location increases with T. 

The effective description of such a transition invokes 
a single coarse grained macroscopic degree of freedom 
~ dislocation deformation characterized by an effective 
mass and a potential energy with two minima. These 
quantities are scaled as some positive powers of L even 
at T 7^ 0, so that the amplitude of the transition between 
the minima decays exponentially as i — >■ cxj - very much 
like d > I systems undergoing first-order transition. 
However, due to the strongly interacting and effectively 
long-range nature of the rough phase, specifics of such 
size dependencies cannot be derived analytically, and we 
have evaluated them numerically. 

The model and its Monte-Carlo simulations. Superclimb- 
ing dislocation is modeled as a quantum string oriented 
along the x-axis and strongly pinned at its both ends 
X = 0,L [l^. The string displacement y{x,t) along the 
y-axis depends on the time t and is measured in units of 
the inter-atomic spacing (« Burger's vector b) with re- 
spect to its equilibrium y = (no tilting is considered). 
The Peierls potential induced by the crystal is taken as 
Up = — Up cos {2TTy{x,t)). The partition function Z has 
the form g, [13 
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+ ^(p-yf + y((V*y)^ + T^,^(V.y)^) 
- iipCOs(27ry(a;,t)) - F?/(x,t)], 



(2) 



where all the variables are periodic in the imaginary time 
t > with the period /3 = 1/T (units h ^ 1, Kb = 
1); the core density p and the superfluid phase are 
canonically conjugate variables, with p' = p — y being the 
local superfluid density; the derivatives Vt^xV, '^t,x4> ^re 
understood as finite differences in the discretized space- 
time lattice (with 200 time slices and x = 0,1,2,...,L 
in units of 5), with Vt,x4> defined modulo 2-k (in order 
to take into account phase-slips); uq^pq stand for the 
average filling factor ( we choose no = 1) and the bare 
superfluid stiffness, respectively, with the bare speed of 



2 



first sound taken as unity. 

Tlie first two terms in Eq.(l2|) describe tlie superfluid 
response of the core and the third term accounts for 
the superclimb effect Q - building the dislocation edge 
so that the core climb y — >■ y ± 1, ±2, ... becomes possible 
by delivering matter p — !■ p± 1, ±2, respectively, along 
the core [H. The dislocation is assumed to be attached 
to large superfluid reservoirs at both ends, with spatially 
periodic boundary conditions for the supercurrent. 

The terms cx m in Eq.Q account for the elastic re- 
sponse of the string, with m and Vd standing for the ef- 
fective mass of the dislocation core (per b) and the bare 
speed of sound, respectively. Since the main source of 
kinetic energy are supercurrents, we have left out the 
term ^ (Vty)^ in Eq.(l2]). The parameter m in Eq.([2]) 
is not actually a constant. It contains a contribution 
from the Coulomb- type interaction potential oc l/\x\ 
between jog s (or kinks, cf. [l5|) separated by a dis- 
tance X [lO|. Accordingly, m has a logarithmic diver- 
gent factor with respect to a wave vector q along the 



core m(q) ~ niQ 



1 + Uc\n(l 



where mp is 



of the order of the atomic ^Hc mass and Uc ~ 1 is a 
parameter characterizing the strength of the interaction 
[1^ [l5| . In solid ^He, the zero-point fluctuation parame- 
ter K = Trh/{AmobVd) ^ 1 [HI. We present our numeri- 
cal results for Uc = 1, V^^ = 5, A' 1. It is important to 
note that the main results are not qualitatively sensitive 
to the long-range interaction. 

The linear force density F ha (ignoring spatial in- 
dices) in Eq.® is determined by the external stress a 
induced by the chemical potential difference Sp, applied 
in the setup @: Sfi fa Spb^, with Sp being the result- 
ing overpressure. Thus a ~ Sp./b'^, i.e. in units of b 
a = S/j, = Sp = F. 

Monte Carlo simulations have been conducted with the 
Worm Algorithm (WA) [l^ for the superfluid part of the 
action, with the Peicrls term treated within the Villain 
approximation similarly to Refs. [l^ . [l5| . The renor- 
malized superfluid stiffness ps (T, F) and compressibility 
k(T,F) have been calculated in terms of the windings 
of the dual variables No significant effect of 

the bias F was found on p^ (T, F) , and thus we consider 
Ps{T,F) = Ps{T,0) = ps{T). In contrast, k{T,F) does 
experience quite dramatic renormalization, which is the 
focus of the present work. 

The bare stiffness po(T) vanishes above some tempera- 
ture Tq comparable to the bulk A-temperature (~ 1 — 2K). 
We have used Tq = 0.2 (in the dimensionless units as 
a fraction of the Debye temperature Tq for the first 
sound), and considered low temperatures - such that 
Ps{T) stayed unchanged within 1-10% of its T = value. 
In other words, the thermal length Lt ~ Ps{0)/T, above 
which Ps{T) becomes suppressed, is the largest scale in 
the problem. 

We have also calculated full xi = SN/5F and dif- 
ferential X2 = dN/dF isochoric compressibilities of the 
dislocation UM. where 5N is the full amount of mat- 
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FIG. 1: (Color online) Renormalized superfluid stiffness Ps{T) 
and the velocity Vs {T, F) of first sound normalized by their 
respective low-T values for different F (shown on the inset), 
L — 30, Up = 3.0. Inset: the region of the dip (cf. Fig. 4 of 
Ref. 01) showing its shifting with F. 



ter accumulated in the extra plane forming the super- 
climbing (edge) dislocation in response to the varia- 
tion of F. These compressibilities obey the relationship 
X2 = Xi + Fdxi/dF, and in the linear regime xi ~ X2- It 
is found that k oc X2/L in the strongly fluctuating regime. 
However, despite such similarity, k and X2/L are not one 
and the same quantity: while X2 describes the shifting of 
the dislocation position, k accounts for the time depen- 
dent response of the superfluid phase (j>. In other words, 
K and Ps enter as renormalized coefficients of the effective 
superfluid action 5^ = J dx J dT[ps{dx4>)'^ /2+k {dt4>)'^ /2] 
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where the renormalized speed of first sound 



(3) 



In the absence of the Peierls potential or at high 
T (where still Ps{T) « Ps(0)), Xi,2 are practically 
equal to the free string response xo ~ L'^/[12mo(l + 
2UcHL/b))V^f] (considered in Rcf. [H for Uc = 0). 
We introduce the normalized quantities Ri^2 = Xi,2/xo- 



i?2 = i?i + FdRi/dF, 



(4) 



so that i?2 = ^1 = 1 for up = 0. 

Our main findings presented below are the following: 
i) a narrow dip in Vs{T, F) vs T at some macroscopically 
small F = Fc, ii) periodicity of Vs{T,F) with respect to 
the bias F; iii) exponential scaling of the dip depth with 
L; iv) hysteresis developing beyond a certain length Lh 
growing with T. 

Dip in the flow rate. The speed ([3]) sets the value of 
the critical flow rate Vcr ~ Vs {T, F) as a threshold for 
generating phase slips in Id - similar to the scenario of 
Ref . [isj . Since k (and R2) exhibits a resonant-type peak, 
Vs acquires the dip shown in FiglTJ Its depth depends 
on L and how closely F is tuned to the threshold value 
Fc Qc L '', 7 PS 1 — 1.7, for the jog-antijog pair creation at 
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FIG. 2: (Color online) Typical behaviors of k, R\^2 vs F: 
L = 56, T /To = 0.05, up = 3.0. Dashed line - the prediction 
of the free string model {up = 0) p^ . 
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FIG. 3: (Color online) Vs{F) normalized as in FiglU The 
parameters are the same as in Fig[2] The narrow dips occur 
at the thresholds F = Fc{L,n), n = 1,2,3 for n jog-antijog 
pairs creation (see the text - Periodicity vs external bias). 



given T (see the inset, Fig. [T]). We believe such a dip has 
actually been detected by Ray & Hallock (cf. Fig. 4 in 
Ref.0). Away from the dip, Vs is practically insensitive 
to F, and the responses k, R2 become essentially linear 
in F [111. 

Periodicity vs external bias. One of the most strik- 
ing features we have found is the quasi-pcriodicity 
in F of K, i?2 as shown in Fig. [2l and in Vs as 
shown in Fig. [31 We attribute this to the reaching 
of the thresholds for creating multiple jog-antijog 
pairs. Thus, the peak (dip) positions are given as 
F = Fc{L,n) sa nFc{L,l) ~ ni-T, n = 1,2,3,..., and 
we predict that the dip in the flow rate observed in 
Ref.[3| (see Fig. 4 there) should recur as a function of the 
applied bias, provided T is kept fixed. 

Size dependencies. Transformation between the smooth 
(where Ri << 1) and the rough {Ri 1) states occurs 
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FIG. 4: (Color online) Height of the first resonance peak 
in R-z(F = Fc) oc K vs L for various parameters (symbols) 
and its fit (lines) by R2 = exp{A + L/Lr) + B with three 
adjustable parameters A,Lr,B. Deviations from the fit at 
large L (shown by arrows) mark the beginning of the hys- 
teretic behavior. Dashed line - the free string model value 
(i?i = i?2 = 1) 
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FIG. 5: (Color onUne) Plots of Fc, Fu, Fl vs L (symbols) and 
their fits (lines). The vertical arrow separates the resonant- 
type peaks in Fig[2] ( L < Lh), and the hysteresis (L > Lh). 
Inset; a typical scan in F shown by the horizontal arrows, 
with the vertical arrows indicating the upper F„ and the lower 
Fl "coercivity" stresses. 



within an exponentially narrow region 6F around Fc. 
It is given by the tunneling rate ~ exp(— L/L/j) << 1 
through the macroscopic jog-antijog barrier (cf. the 
mechanism for the kink-antikink tunneling, Ref . [TTI| ) . 
where Lr stands for the tunneling length. Thus, as fol- 
lows from Eq.(|3]), the peak value of R2{T,Fc) diverges 
with L as w l/SF - exp(i/L/?). Fitting R2iT,Fc) by 
an exponential function, FigUl gives « — 
T/Tnf-^. Here is the T = tunneling length 
[Lq « 1.7 for Up = 3.0) and Tr sets the scale for ther- 
mal roughening, i.e. the temperature above which the 
density of jog pairs is large even in the limit _F — > 0. 
Tji is determined by the double energy 2A oc y/up of a 
jog modeled as a Sine-Gordon soliton. We have found 



4 



Tr oc Up, s = 0.5 ± 0.1 which is consistent with such an 
interpretation. 

The critical stress Fc{L, 1) has been found to deviate 
from the l/L dependence [11| at finite T. Specifically, 
Fc cx with 7(T) > 1. As T grows, -/{T) 1.7, 

and 7(r) — > 1 in the limit T ^ 0. The temperature 
scale for this variation is set by the Peierls potential 
amplitude up as well. It is natural to attribute the 
deviation from 7 = 1 to a suppression of the energy gap 
A oc L^-'^ at finite T. 

Hysteresis. The resonant-peak type behavior in i?2 (and 
in k) turns out to be a precursor for the jump in Ri. As 
seen from Eq.g]), R2 « FcdRi/dF - Fc'/SF as SF 0. 
The hysteresis emerges when SF becomes significantly 
less than the coexistence region for the smooth and rough 
states. This conditions sets a typical length Lh above 
which {L > Lh) hysteresis develops. We have found that 
Lh grows with T as Lh ^ Lq{T /Tg)'^'^ ^ Lo, 7/i > 
(7/1 « 2 — 3 for Up = 1 — 3). We attribute the en- 
ergy scale Tg to the tunneling splitting energy through 
the microscopic jog-antijog barrier so that Tg ^ Tr ^M- 
This feature is clearly due to the collective multi-jog na- 
ture of the rough state, and it deviates strongly from the 
single pair tunneling scenario [ll| (see Eq.(39)) where 
the tunneling rate saturates at some length decreasing as 

FigOdemonstrates the L-dependencies of the "coerciv- 
ity fields" .Along the lower branch of the hysteresis loop 
(inset in Fig|5]) the dislocation is in the smooth state. 
Upon increasing _F, it "jumps" into the rough state at 
F Fu- While "moving" back along the upper branch 
representing the rough state, the dislocation returns into 
the smooth state aX F k. Fl < Fu- As seen from the main 
panel of FigO both fields scale as F^^l oc L"'''"'^, with 
7„ « 0.4 and 71, « 2.7, respectively. The middle straight 
line corresponds to the extrapolation of the peak posi- 
tion data Fc oc L~'^ , n = 1, 7 ~ 1.7 (for L < Lh). Such 



strong sensitivity to the size L as well as Lh growing with 
T clearly indicate that the stress-induced roughening is 
a phase transition at finite T in Id. 

The resonant-type behavior and the hysteresis have 
also been found in simulations of gliding dislocation 
p^ . While Lt determines the upper spatial scale for 
the supcrclimb, no such restriction exists for the glide, 
so that the formal limit L ^ 00 can be considered. 

Discussion and conclusions. The narrow dip in the su- 
perfiow rate observed in Ref. Q may have its origin in the 
stress induced roughening effect of superclimbing dislo- 
cations. Specifically, biasing a superfiuid dislocation net- 
work by macro scopically small overpressure can induce 
strong suppression of the first sound along the superfiuid 
cores. Such suppression is characterized by the periodic- 
ity of the dip in the fiow rate, FiglSl which can be used 
as the experimentum crucis for the proposed scenario. 
We estimate the critical overpressure Sp in terms of L 
and a typical jog-antijog energy 2 A « O.IK (cf. [Tsj ) 
as Sp/p w 2Ab/{LTD) ^ lQ~'^b/L. Measurements of Sp 
where the dip recurs can provide crucial information on 
the nature of the dislocation network - its typical free 
segment length L. 

As T is lowered, the opposite condition L > Lh 
is fulfilled so that the dislocation behavior becomes 
hysteretic between its smooth and rough states. We 
believe it is also important to study the hysteresis in the 
fiow rate vs chemical potential (the upper F„ and the 
lower Fl fields) at different temperatures. 
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